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Abstract 

Magnetic properties of single crystals of the heterometallic complex 
[Cu(bpy) 3 ]2[Cr(C 2 04)3]N03-9H 2 (bpy = 2,2'-bipyridine) have been 
investigated. From the recorded EPR spectra, the spin-Hamiltonian 
parameters have been determined. The magnetization measurements 
have shown magnetic anisotropy at low temperatures, which has been 
analysed as a result of the zero-field splitting of the Cr m ion. By 
fitting the exactly derived magnetization expression to the measured 
magnetization data, the axial zero- field splitting parameter, D, has 
been calculated. Comparing to the EPR measurements, it has been 
confirmed that D can be determined from the measurements of the 
macroscopic magnetization on the single crystals. 

1 Introduction 

In the past several years a great number of homo- and heteropolynuclear 
complexes containing different transition metal ions have been synthesized 
and their magneto-structural properties have been investigated. Interaction 
between paramagnetic centres involved as well as the crystal structures of 
these compounds determine their magnetic properties. In describing mag- 
netic behaviour of these species, besides the exchange interaction between 
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metal ions, it is often necessary to include the ligand field effects which induce 
the anisotropy of paramagnetic centres and zero-field splitting. Paramagnetic 
metal centres of these complexes are usually bridged by diamagnetic ligands 
which mediate exchange interaction between metal ions pp. In that sense, 
a widely used ligand is the oxalate C^O 2- anion - not only because of the 
influence it has in transmitting electronic effects between magnetic centres, 
but also because of its extraordinary rich binding facilities [2]. The use of 
stable mononuclear anionic oxalate complexes, [M m (C 2 04)3] 3_ (M m = Cr, 
Fe, Ru), as ligands toward other metal ions provides an efficient route for the 
synthesis of heteropolynuclear oxalate-bridged species of different nuclearity 
and dimensionality, and with a variety of magnetic properties [3-5]. 

In this paper magnetic properties of the heterometallic complex [Cu(bpy) 3 ] 2 - 
[Cr(C 2 04)3]N03-9H 2 (bpy = 2,2'-bipyridine) in the form of single crystals 
have been investigated. The compound, in which the oxalate anion does not 
act as a bridging ligand but as a terminal one, was synthesized by the reac- 
tion of the [Cr(C 2 4 )3] 3 ~ and [Cu(bpy) 3 ] 2+ complex ions as building blocks, 
during the course of a more extensive research work on oxalate-based transi- 
tion metal species [B] . The magnetic measurements performed previously on a 
powder sample of this heterometallic complex showed a small deviation of the 
effective magnetic moment from the spin-only value for the three uncoupled 
spins (Sea, Sqt, Scu) = (1/2, 3/2, 1/2) at low temperatures. This observation 
could not be well understood in the model of exchange interaction between 
Cr 111 and Cu ions. It was suggested that the low temperature magnetic 
properties of this complex could be explained by considering the zero-field 
splitting on Cr 111 . For this purpose, in the present study, the magnetic mo- 
ment of the single crystals of the complex has been measured. Also the EPR 
measurements have been performed in order to obtain the spin-Hamiltonian 
parameters. Magnetization measurements have been analysed with respect 
to the crystal structure of the compound and the EPR measurements using 
an exact, instead of the Van Vleck, approach. 

2 Experimental 

Electron paramagnetic resonance (EPR) and magnetization measurements 
were performed on the single crystals and also on a powdered sample of 
[Cu(bpy)3] 2 [Cr(C 2 4 )3]N03-9H 2 0. Transparent-blue single crystals of the 
complex were obtained by slow evaporation of an aqueous solution contain- 
ing the [Cr(C 2 4 )3] 3 " and [Cu(bpy) 3 ] 2+ ions in the ratio 1:1 [BJ. Dimensions 
of the prepared plate-like single crystals were approximately 2.0 x 0.4 x 0.2 
mm 3 up to 3.0 x 0.5 x 0.5 mm 3 . 
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For the EPR measurements, the crystal was mounted on a quartz holder 
in the cavity of an X-band EPR spectrometer (Bruker Elexsys 580 FT/CW) 
equipped with a standard Oxford Instruments model DTC2 temperature 
controller. The measurements were performed at the microwave frequency 
around 9.6 GHz with the magnetic field modulation amplitude of 3 G at 
100 kHz. The EPR spectra were measured at different temperatures, from 
the room to the liquid helium temperature. It was found that the spectra 
were approximately temperature independent (except for the paramagnetic 
behaviour of the spectral line); therefore, only the room temperature spectra 
are presented in the paper. In the first set of the EPR measurements the 
crystal was mounted to the quartz holder so that the crystallographic b axis 
was perpendicular to the magnetic field. The crystal was rotated around the 
b axis and the EPR spectra were recorded every 5°. Since the complex has 
a monoclinic crystal symmetry, the a* and c* axes, mutually perpendicular 
and perpendicular to the b axis, were chosen arbitrarily. The first rotation 
was done around the b axis in the a* — c* plane. For the second set of 
measurements, the crystal was rotated for 90° in such a way that it performed 
rotation in the a* — b plane and the third set included the rotation in the 
b — c* plane. Hence, all three rotations were performed around the chosen 
laboratory axes: a*, b and c* with the magnetic field in the plane of rotation. 
The rotation was controlled by a goniometer with the accuracy of 0.5°. The 
largest uncertainty (2-3°) was related to the optimal deposition of the crystal 
on the quartz holder. 

For the magnetization measurements, a single crystal of the complex was 
attached to a homogeneous straw using a small dot of vacuum grease. Mag- 
netic moment of the sample was measured using the commercial SQUID 
magnetometer MPMS-5 (Quantum Design). Measurements were performed 
in the temperature range 1.8-300 K for the applied magnetic field of B = 1 
T. Also, magnetic moment of the sample in the dependence on the applied 
magnetic field 0-5.5 T at several temperatures was measured. The measure- 
ments were performed for two different orientations of the single crystal with 
respect to the applied magnetic field: for the magnetic field parallel to the 
crystallographic b axis (B \\ b) and for the magnetic field perpendicular to the 
crystallographic b axis (B _L b). Two crystals were used during the measure- 
ments. Mass of the first crystal (for B \\ b) amounted 0.136 ± 0.001 mg and 
mass of the second crystal (for B _L b) amounted 0.070 ± 0.005 mg. Because 
of the very small masses of the crystals, measurements of the temperature 
dependence of the magnetic moment were carried out in the magnetic field of 
1 T in order to obtain a well measurable signal. The calculated molar magne- 
tizations were corrected for the diamagnetic contribution of the constituent 
atoms, temperature independent paramagnetism of the Cr m ion and two 
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Cu 11 ions [7] and for a tiny contribution of the background (vacuum grease). 

3 Results and Discussion 
3.1 Structural considerations 

The details of crystal structure of the complex [Cu(bpy) 3 ] 2 [Cr(C204)3]N03 — 
•9H 2 (bpy = 2,2'-bipyridine) were described previously [Bj. The complex 
crystallizes in the monoclinic space group P2x/c, with the unit cell parame- 
ters a = 31.314(3) A, b = 13.5356(11) A, c = 22.202(2) A, (3 = 132.012(13)° 
and Z — 4. The paramagnetic centres (Cr m in the [Cr(C204)3] 3_ anion and 
two Cu 11 in two [Cu(bpy) 3 ] 2+ cations) are located in distorted octahedral 
coordination environments. Distortions are expected to be trigonal (due to 
the rigidity of the didendate oxalate and bypiridine ligands, respectively) 
and additionally, in the case of [Cu(bpy)3] 2+ cations, tetragonal (due to the 
Jahn- Teller effect). The tetragonal component of the [Cu(bpy)s] 2+ cation dis- 
tortion can be easily recognized in the differences of the Cu-N bond lengths 
as described in [6J . For the ligand field environment of the Cr m ion the tetrag- 
onal distortion is negligible as seen from the very small fluctuations of the 
Cr-0 bond lengths (~ 0.02 A). Instead, rather equal values of the three bite 
angles (~ 82.7°) show a trigonal type of distortion; bite angles are defined 
as the angles between two coordinated oxygen atoms from the same oxalate 
ligand and the Cr m ion. Geometrical parameters for the trigonal distortion 
(compression ratio, s/h and angle of twisting, cf) [S]) can be estimated from 
the positions of six coordinated oxygen atoms: 01, 02, 05, 06, 09 and O10 
(Fig. 1(a)). Assuming that the [Cr(C204)s] 3 ~ anion has an approximate D% 
symmetry and that a local three-fold rotation axis passes through the centres 
of the triangles (01, 05, 09) and (02, 06, O10), the values for s/h and are 
found to be 1.32(3) and 52.2(7)°, respectively, quite different from the ideal 
octahedral values of 1.22 and 60° [8]. The mentioned three-fold rotation axis 
is the only one rotation symmetry of the [Cr(C204)3] 3_ anion which transfers 
the oxalate groups into each other. 

As a consequence of the symmetry operations in the P1\jc space group, 
four [Cr(C 2 4 )3] 3 ~ anions are located in one unit cell (labeled as I, II, III 
and IV in Fig. 1(b)). Local three-fold rotation axes of [Cr(C204)3] 3_ anions 
are inclined for 9 = 11.77(1)° with respect to the crystallographic b axis. The 
anions I and IV (as well as the anions II and III) have identical orientation in 
the space because they are related through the centre of inversion i (shown 
in Fig. 1(b)). The anions related through the two- fold screw symmetries 
(I and III or II and IV) have mutually different orientation in space as a 
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(a) 




(b) 

Figure 1: (a) The local, approximate three-fold rotation axis of the 
[Cr(C 2 4 )3] 3_ anion (white rod with a triangle on the top) which passes 
through the Cr m ion and through the centres of two triangles of coordinated 
oxygen atoms: (01, 05, 09) and (02, 06, 010). (b) Four symmetry related 
[Cr(C204)s] 3 ~ anions (I, II, III and IV) in the unit cell of the investigated 
compound showing the inclination angle 9 of their local three-fold rotation 
axes with respect to the crystallographic h axis. Two two-fold screw sym- 
metry axes and the centre of inversion % are also shown. The [Cu(bpy)s] 2+ 
cations and other constituents of the complex are omitted for the reasons of 
clarity. 

consequence of the above mentioned angle of 9 = 11.77(1)° between the local 
three-fold rotation axes and the crystallographic b axis (b axis is parallel with 
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the two- fold screw symmetry axes). Therefore, the local three-fold rotation 
axes of the anions I and IV make an angle of 23.54(1)° with respect to the 
local three-fold rotation axes of the anions II and III (Fig. 1(b)). 

3.2 EPR study 

The EPR spectra of a single crystal of the [Cu(bpy) 3 ]2[Cr(C204)3]N0 3 -9H 2 
complex when the magnetic field was parallel to the a*, b and c* axes are 
shown in Fig. El In the spectra two types of signals can be separated: the 
first one from Cu 11 appearing always near 3232 G and the second one from 
Cr m which shows a large angular dependence: from ~ 1700 — 1800 G (B \\ a* 
and B || c*) to ~ 3300 G (B || b) and even splitting into two lines when the 
crystal was rotated around the a* and c* axes. The third signal, at 1750 G, 
comes from the EPR cavity. 
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Figure 2: The EPR spectra of a single crystal of the investigated complex 
recorded at room temperature, with the magnetic field B parallel to the 
laboratory a*, b and c* axes. 

Fig. [3 presents the experimental and computer simulated angular depen- 
dence of the resonant field for the single crystal of the complex for all three 
rotations, i.e. for the rotations around the c* (from B \\ b to B || a*), b (from 
B || a* to B || c*) and a* (from B \\ c* to B \\ b) axes. A computer simulation 
was carried out by the XSophe-Sophe-XeprView software [9]. 
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Figure 3: Experimental (circles and squares) and simulated (solid lines) an- 
gular variation of the resonant field for the single crystal of the investigated 
complex at room temperature. 

The observed spectral lines are very broad (~ 350 G) mainly because of 
the dipole-dipole interaction so the spectra from Cu 11 have been explained 
approximately by the Zeeman interaction of the electron spin S — 1/2 with 
the magnetic field: 



g-tensor and S is the electron spin operator. It has been assumed that the 
g-tensor for Cu 11 is isotropic with the ^-factor value g = 2.11 ± 0.02. In 
this paper, we have not further analysed the anisotropy of the g-tensor for 
Cu 11 because its resonant field, compared to that from Cr m , depends on 
orientation only weakly. 

For Cr m , the following form of the spin-Hamiltonian has been assumed: 



where D is the zero-field splitting tensor. The relation (2) could be written 
equivalently (for more details, see for example [10]): 




H = pi B B • g • S + S • D • S 



(2) 



H = hbB-s-S + D S z 2 - 



S(S+1) 
3 



+ E (S x 2 - S y 2 ) 



(3) 
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where D and E are the axial and rhombic zero-field splitting parameters, re- 
spectively, and the other parameters have their usual meanings. We have 
assumed an isotropic g-tensor for Cr m . A set of parameters reproduc- 
ing experimental data the best is the following one: g = 1.963 ± 0.002, 
D = (0.63 ±0.01) cm" 1 , \E\ = (0.02 ±0.01) cm" 1 and the angle between the 
magnetic z axis and the laboratory z axis (which corresponds to the crystal- 
lographic b axis) is 9 = (11.5 ± 0.5)°. This angle is in excellent agreement 
with the previously mentioned value of 11.77(1)° for the angle between the 
calculated three-fold rotation axes of Cr m and the crystallographic b axis 
(Sect. 3.1). Two signals from Cr IH observed in the b — a* and c* — b planes 
(Fig. 3) occur because two orientations of the three-fold rotation axes ex- 
ist (Fig. 1(b)). The EPR transition for Cr m which can be seen in Fig. |2] 
corresponds to the transition Ms = — 1/2 <-> Ms = 1/2. Another allowed 
transition Ms = —1/2 «-> Ms = —3/2 should be observed at a higher mag- 
netic field which was out of range of the experimental set-up (behind 1 T) and 
the transition M s = 1/2 <-> M s = 3/2 could not be observed by the X-band 
EPR spectrometer, because the zero field splitting D = 0.63 cm -1 pa 20 GHz 
is much higher than our microwave frequency of 9.6 GHz [TU]. These EPR 
measurements could give only an absolute value of the parameter D, but the 
positive sign was chosen in accordance to the magnetization measurements. 
The sign of E has no physical meaning, except in terms of convention of the 
chosen a* and c* axes [11]. 

In order to additionally check the spin-Hamiltonian parameters obtained 
from the single crystal data, the EPR spectrum of a powdered sample of the 
complex is recorded at room temperature, Fig. HI 

The powder simulation is performed using the previously mentioned pa- 
rameters for Cu 11 and for Cr m obtained from the single crystal simulation. 
The total simulation line of the powder spectra is obtained as a sum of the 
simulation lines for Cu 11 and for Cr 111 in a ratio 2 : 1 according to the fact 
that there is two times more Cu 11 than Cr m ions in the unit cell. It could be 
seen from Fig. H] that the simulation is in good agreement with the powder 
spectrum. 

3.3 Magnetization study 

The results of magnetization measurements of the [Cu(bpy)3]2[Cr(C204)3]- 
N03-9H 2 complex for the parallel and perpendicular orientations of the 
crystallographic b axis with respect to the applied magnetic field B = 1 T, 
i.e. for B \\ b and B ± b, are shown in Figs. and El respectively, as the 
dependence of the product of molar magnetization and temperature, M mo {T 
on temperature T. 



8 




experiment 

simulation Cu 1 



simulation Cr 

total simulation 



1000 



2000 



3000 



4000 



5000 



B[G] 



Figure 4: The powder EPR spectrum of the investigated complex recorded 
at room temperature (thick solid line). The simulations for Cu 11 (dashed 
line) and for Cr 111 (dotted line) are obtained using the single crystal param- 
eters. The total simulation line (thin solid line) is obtained as a sum of the 
simulation lines for Cu 11 and for Cr 111 in a ratio 2:1. 

At high temperatures the value of the product of molar magnetization 
and temperature, M mo {F should correspond to the spin-only value. For the 
three uncoupled spins («Scu, So-, Sen) = (1/2, 3/2, 1/2) this value is given by 
the expression: 



M mol T = -^pL[g 2 Cu Scn(S C u + 1) + g 2 Cl S C r{Scr + 1) + <£ u S C u(S C u + 1)] (4) 



where Na is Avogadro's number and k is the Boltzmann constant. By using 
the ^-factors gc r = 1.963 and gc u = 2.11 (obtained from the EPR measure- 
ments) the M mo iT product equals 2.642 • 10 4 emu mol -1 K (shown as the 
dashed line in Figs. [5] and [6]). The value of the M mo {T product for B || b 
(Fig. E]) is constant at high temperatures and corresponds to the predicted 
spin-only value. The M mo iT product increases on cooling, attains a maxi- 
mum value at 5 K and thereupon decreases. For B _L b (Fig. [6]), the M mo {T 
product is also constant at high temperatures, having the value as expected 
for the three independent spins. On cooling, this value remains constant 
down to 15 K and decreases on further cooling. The relative deviation of the 
M mo {T value at low temperatures from the spin-only value is rather small, 
i.e. approximately 9% for B || b and 17% for B _L b. 
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Figure 5: The M mo {F versus T plot for the parallel orientation of the crys- 
tallographic b axis with respect to the applied magnetic field B = 1 T for 
the investigated complex. The horizontal, dashed line corresponds to the 
theoretical value expected for the three uncoupled spins. 

In the previous investigation of magnetic properties of this heterometal- 
lic complex [6], it was observed that the long-range intermolecular interac- 
tion was negligible. Also, exchange interaction between Cu 11 and Cr m ions 
and between two Cu 11 ions can be discarded because in the investigated 
heterometallic complex Cu 11 and Cr m ions are not bridged by diamagnetic 
ligands which could mediate the exchange interaction. Besides, by using 
the model of exchange interaction between paramagnetic centres it was not 
possible to explain the decrease of the M mo {F product at low temperatures 
and the observed magnetic anisotropy [6J. Therefore, in further analysis of 
magnetic properties of the complex we have excluded exchange interaction 
between the magnetic ions. Considering the EPR measurements, we have 
tried to associate the observed magnetic anisotropy at low temperatures to 
the zero-field splitting of Cr m ion levels in a distorted octahedral surrounding 
((^-factors for Cr m and Cu 11 ions assumed to be isotropic). 

With this aim in view, in this paper we have derived the exact expression 
for the molar magnetization of the complex. The Van Vleck approach has 
not been used because the condition jisB/kT <C 1 is not fulfilled at low 
temperatures and for the applied magnetic field 1 T. Moreover, the deviation 
of the M mo {T product from the spin-only value appears at relatively low 
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Figure 6: The M mo \T versus T plot for the perpendicular orientation of the 
crystallographic b axis with respect to the applied magnetic field B = 1 T 
for the investigated complex. The horizontal, dashed line corresponds to the 
theoretical value expected for the three uncoupled spins. 



temperatures and is relatively small. Also, with regard to the EPR results 
and other studies of similar heterometallic complexes in which Cr m ions also 
adopt distorted octahedral coordination [T2J, [13], it is reasonable to expect 
that the value of the axial zero-field splitting parameter is of the order 1 cm -1 . 
Thus, neither the condition D ^> gcr^BBj_, usually used in the Van Vleck 
formula, is valid for the applied magnetic field of 1 T. In literature where 
zero-field splitting is studied, commonly the Van Vleck formula is used with 
an additional energy linearization pQ. In the present investigation the idea 
has been to find out to what extent the approximations used in the Van Vleck 
approach have the influence on the interpretation of magnetic measurements 
and on the quantitative determination of parameter D, too. 

Generally, molar magnetization of a quantum system with an energy spec- 
trum E n , n = 1,2,3,... in the presence of a magnetic field B in thermal 
equilibrium is given by the Boltzman distribution law: 

AT, V dE n -En/kT 
jv ^Zj dB e 

Mm ° l = £ e -E n /kT ( 5 ) 

n 

The equation (jHJ) is the fundamental expression for magnetization and does 
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not lean on any approximation. Energy spectrum E n , n = 1, 2, 3, . . . corre- 
sponds to the eigenvalues of the Hamiltonian of the given system. 

The Hamiltonian of the Cr m ion in a distorted octahedral surrounding in 
the presence of the magnetic field is given by the relation (j3J), which can be 
written in an equivalent form: 



H = gcv^B S z • Bn + gcrf^B S x -B± + D 



S(S + 1) 



+ E(S y 



S 2 



) (6) 



where B\\ and B± are components of the applied magnetic field parallel and 
perpendicular to the anisotropy axis of the local surrounding of the Cr m ion, 
respectively, the g-tensor is assumed to be isotropic with the go principal 
value and the other symbols have their usual meanings. The energy lev- 
els E n cr are determined by diagonalization of the above Hamiltonian. The 
components of molar magnetization parallel My and perpendicular M± to 
the anisotropy axis are given by the following equation: 



V 9E n,Cv -E nC .r/kT 
^ dBn I 



M,u ;Cr = JVa— 



J2 e-^.cr/fcT 



(7) 



Generally, when the magnetic field B is applied in an arbitrary direction 
(i.e. neither parallel nor perpendicular to the anisotropy axis), the resultant 
magnetization M is not parallel to the applied field B. The magnetic field 
vector B and the magnetization vector M with its parallel and perpendicular 
components for such a general case are shown in Fig. [7J 

The total magnetization is a vector sum of the Mm and Mj_ components. 
As presented in Fig. [7J 5 is the angle between the applied magnetic field B 
and the anisotropy axis z and ip is the angle between the total magnetization 
M and the anisotropy axis z. Projection of the total magnetization to the 
direction of the applied magnetic field is measured in the experiment (M meas ). 
Therefore, contribution of the Cr m ion to the molar magnetization of the 
complex is given by the following expression: 



M Cr (B,T,6) 



1VI \\,Cr ^ 1VI ±,Cr 



COS 



5 — arctan 



M 



_L,Cr 



Mm 



Cr 



(8) 



In addition to the contribution of the Cr m ion to the molar magnetization 
of the compound, there is also a paramagnetic contribution of the two Cu 11 
ions. The Hamiltonian of the Cu 11 ion in an external magnetic field is given 
by (00) and its eigenvalues E ni c u are used to evaluate molar magnetization of 
the Cu 11 ion according to the equation (jSJ). 
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Figure 7: Schematic view of the magnetic field and magnetization vectors, 
B and M, respectively, when B is applied in an arbitrary direction to the 
anisotropy axis z. B\\ and B± are components of B (Mm and M± are compo- 
nents of M) parallel and perpendicular to the anisotropy axis, respectively; 
M is the resultant magnetization and M meas is a projection of the magneti- 
zation to the direction of the applied field measured in the experiment. 

In the experiment the magnetic field B was parallel or perpendicular to 
the crystallographic b axis. In these two specific cases both anisotropy axes 
(i.e. local three-fold rotation axes) form the same angle 5 with the magnetic 
field, as it can be observed in Fig. Q](b). The final expression for the measured 
magnetization when the magnetic field is parallel or perpendicular to the 
crystallographic b axis is: 

M meas (B } T } 6) = M C r(B,T,5) + 2M Cu (B 1 T). (9) 

The nonlinear fitting of the measured M mo {F data versus temperature de- 
pendence in the temperature range 1.8-60 K for both crystals was performed 
by the equation ([9]) multiplied by temperature T. Fitting parameters were 
the axial zero-field splitting parameter D, the angle between the anisotropy 
axes and the magnetic field 8 and a multiplication constant taking into ac- 
count uncertainty of the weighed mass of the sample. The values of the 
g-factors were fixed to go = 1-963 and gcu = 2.11 and the rhombic zero-field 
splitting parameter to E — 0. The nonlinear fitting gave the following values: 
D = 1.03 cm -1 , 5 = 69.9° for the parallel orientation and D = 0.97 cm -1 , 
5 = 48.2° for the perpendicular orientation. The discrepancy factor, defined 
as R = £ [(M mol T) obs - (M mol T) calc ] 2 /Z[(M mo iT) obs ] 2 equaled 5.7 • lO" 6 
and 9.0 • 10~ 6 for the parallel and perpendicular orientations, respectively. 
The obtained parameters D are consistent for both orientations. According 
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to the definiton and geometrical reasons, the angle 5 in the B _L b case should 
be in the interval [90° — 5(B || b), 90°]. The obtained values of the angle 5 
in both cases fulfill this requirement. The calculated curves for B \\ b and 
B _L 6 are presented in Figs. [5] and El respectively, showing good overlapping 
with the measured data. 

Several additional measurements of the temperature dependence of the 
magnetic moment of the complex for B \\ b and B _L b were performed on 
different single crystals. The nonlinear fitting gave mutually consistent values 
of the parameters D and S. 

For the B || b case the angle 5 should correspond to the previously men- 
tioned angle between crystallographic b axis and the anisotropy axis 9, which 
is determined from the crystallographic and EPR data (Section 3.1 and 3.2). 
The calculated curve with the parameters obtained by the nonlinear fitting 
of the M mo iT data in the B \\ b case with the value of the angle 5 fixed to 
the value of 9 = 11.77° was not in accordance with the measured data. So, 
the angle 6 is used as one of the fitting parameters. In the B _L b case the 
value of the angle 5 could not be predicted from the structural considera- 
tions. Therefore, in the nonlinear fitting in the B _L b case 5 is also taken as 
fitting parameter. 

The nonlinear fitting of the molar magnetization M mo i versus magnetic 
field B dependence in the range 0-5.5 T was performed for the measurements 
at 2 K. It can be shown that the shape of the M mo i(B) curve is rather 
insensitive to the values of D and 5. So, all parameters as before were used 
except the multiplicative constant. The results of the fit have shown that the 
measured data can be well described by the relation ([9]). The experimental 
data and the calculated curves of M mo i(B) for B \\ b and B _L b are given in 
Fig. [HI They clearly show that a linear approach in magnetic susceptibility 
calculation would be meaningless for the used ranges of magnetic field and 
temperature. 

In the Van Vleck approach (described in [TJ) it is assumed that the en- 
ergy levels E n can be expanded according to the increasing powers of B 
(E n = E { n ] + E { n ] B + E ( n ] B 2 + ■■■) and that the condition B/kT < 1 is 
fulfilled. Introducing these ap The nonlinear fitting of the formula obtained 
by the Van Vleck approach to the experimental data gave about two times 
higher values of the axial zero-field splitting parameter D, which is also in 
bigger discrepancy with the results of the EPR measurements. Fit to the 
experimental data with the parameter E as fifth fit parameter did not give 
physical and reliable results because of overparametrization and other nu- 
merical problems. 

The values of the parameter D obtained from two independent measure- 
ments, that is for two different orientations of the single crystal in the mag- 
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Figure 8: The M mo i versus B plot for the parallel and perpendicular orien- 
tations of the single crystals of the investigated complex with respect to the 
applied magnetic field at 2 K. 

netic field are consistent to each other, and temperature dependence of the 
magnetization of the complex can be well described by this model. Uncer- 
tainty of the obtained parameters can be estimated in the following manner. 
First, it should be emphasized that the effect of zero-field splitting is observ- 
able only at low temperatures and that the deviation of the measured magne- 
tization from the spin-only value is rather small. Therefore, the uncertainty 
of the mass of the sample as well as the inability of the exact determination of 
the contribution of the background to the measured magnetization contribute 
nonnegligibly to the uncertainty of the obtained parameters. The influence 
of the background has been estimated in such a manner that the molar mag- 
netization at high temperatures corresponds to the spin-only value. From 
the nonlinear fitting of the measured data corrected for the slightly different 
values of the background magnetization it has been observed that the value 
of the parameter D is changing within ±0.2 cm -1 and the value of 5 within 
±1°, whereas the discrepancy factor R does not change significantly. 

The EPR measurements confirm the results of the measurements per- 
formed on SQUID magnetometer taking into account associated uncertainties 
of the calculated parameters and the specificities of the used experimental 
methods. EPR measurements are suitable for accurate determination of the 
spin-Hamiltonian parameters of a system. Magnetization measurements can 
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provide only an effective value of the zero-field splitting parameter. Although 
the zero-field splitting of the energy levels of the Cr m ion determines mainly 
the low temperature behaviour of the magnetization of the complex, it is pos- 
sible that there are some other, much smaller effects which also contribute to 
the magnetization and may influence the calculated values of the parameter 
D and angle 5. Hence, the anisotropy of the g-factors of the Cr m and Cu 11 
ions, the exchange interaction between Cr m and Cu 11 ions and between two 
Cu 11 ions are likely present to a certain extent and possibly have an effect 
on the calculated value of the parameters. Besides, the magnetization has 
been measured in the relatively high field (1 T) which may influence the 
magnitude of anisotropy and the direction of the anisotropy axis. Magne- 
tostrictive effects in such a structure are not known, but it could be possible 
that the applied field changes the direction and the amount of anisotropy or 
even produces some additional contribution [TJ]. This could be the reason 
of differing between the magnetization and EPR determinations of the zero- 
field splitting parameter D and of the angle between crystallographic b axis 
and the anisotropy axis 9. 



4 Conclusion 

Magnetic properties of the heterometallic complex [Cu(bpy)3] 2 [Cr(C 2 04) 3 ]- 
N03-9H 2 (bpy = 2,2'-bipyridine) in the form of the single crystal have been 
investigated. At low temperatures magnetic anisotropy has been observed. 
Taking into consideration the results [6] of the previous magnetic measure- 
ments and known crystal structure of the complex it has been concluded 
that the exchange interaction between Cr m and Cu 11 ions and between two 
Cu 11 ions, as well as the long-range intermolecular interaction are negligible. 
The zero-field splitting of the energy levels of the Cr m ion in the trigonally 
distorted octahedral surrounding leads to the magnetic anisotropy at low 
temperatures. This effect also completely determines the low temperature 
behaviour of the magnetization of the complex because other possible inter- 
actions are negligible, unlike to similar heterometallic complexes where other 
interactions, such as exchange interaction between paramagnetic centres or 
long-range interaction dominate. Therefore, this complex is very suitable for 
qualitative and quantitative studies of zero-field splitting. The theoretical 
expression for molar magnetization of the complex has been derived without 
use of any approximation. The nonlinear fitting to the experimental data has 
shown that the magnetic properties of the compound can be well described 
by this model. 

The calculated value of the axial zero-field splitting parameter of the 
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Cr ion is in agreement with those previously reported for similar polynu- 
clear complexes [T2l [TB"] . Also, the sign of the parameter D > has been 
uniquely determined. Besides the magnetic moment measurement of the sin- 
gle crystals of the complex, the behaviour of the EPR spectra by rotating a 
single crystal about the three independent axes has been studied. From the 
spectral line shifting the spin-Hamiltonian parameters have been obtained, 
which are in agreement with the data from the literature [15-17]. It has 
been also shown that the Van Vleck approach gives bigger disagreement be- 
tween two methods. The EPR measurements, which are more relevant for 
microscopic parameters determination, thus have confirmed the parameters 
obtained from the magnetic bulk measurements. It has been shown that 
magnetic measurements are also appropriate for the anisotropy parameter 
determination if appropriately analysed, i.e. if using exact calculation of the 
magnetization of a system. 
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